Fill-in reduction in sparse matrix factorizations using hypergraphs
نویسندگان
چکیده
We discuss partitioning methods using hypergraphs to produce fill-reducing orderings of sparse matrices for Cholesky, LU and QR factorizations. For the Cholesky factorization, we investigate a recent result on pattern-wise decomposition of sparse matrices, generalize the result, and develop algorithmic tools to obtain more effective ordering methods. The generalized results help us to develop fill-reducing orderings for LU factorization in a similar way to those for Cholesky factorization, without symmetrizing the given matrix A as |A| + |A | or |A ||A|. For the QR factorization, we adopt a recently proposed technique to use hypergraph models in a fairly standard manner. The method again does not form the possibly much denser matrix |A ||A|. We also discuss alternatives for LU and QR factorization where the symmetrized matrix can be used. We provide comparisons with the most common alternatives in all three cases.
منابع مشابه
Riordan group approaches in matrix factorizations
In this paper, we consider an arbitrary binary polynomial sequence {A_n} and then give a lower triangular matrix representation of this sequence. As main result, we obtain a factorization of the innite generalized Pascal matrix in terms of this new matrix, using a Riordan group approach. Further some interesting results and applications are derived.
متن کاملAlgorithms for Approximate Subtropical Matrix Factorization
Matrix factorization methods are important tools in data mining and analysis. They can be used for many tasks, ranging from dimensionality reduction to visualization. In this paper we concentrate on the use of matrix factorizations for finding patterns from the data. Rather than using the standard algebra – and the summation of the rank-1 components to build the approximation of the original ma...
متن کاملRandomized Sparse Direct Solvers
We propose randomized direct solvers for large sparse linear systems, which integrate randomization into rank structured multifrontal methods. The use of randomization highly simplifies various essential steps in structured solutions, where fast operations on skinny matrix-vector products replace traditional complex ones on dense or structured matrices. The new methods thus significantly enhanc...
متن کاملDirected Moore Hypergraphs
For graphs with maximum degree A and diameter D, an upper bound on the number of vertices is 1 + AxF:-,‘(A 1)‘. This bound is called the Moore bound for graphs and the graphs that attain it are called Moore graphs. Similar bounds for directed graphs and for hypergraphs have been defined and the existence of directed Moore graphs and of Moore hypergraphs has been studied. In this article, we def...
متن کاملPreconditioned Conjugate Gradient Methods for Three Dimensional Linear Elasticity
Finite element modelling of three dimensional elasticity problems gives rise to large sparse matrices. Various preconditioning methods are developed for use in preconditioned conjugate gradient iterative solution techniques. Incomplete factorizations based on levels of fill, drop tolerance, and a two level hierarchical basis are developed. Various techniques for ensuring that the incomplete fac...
متن کامل